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Abstract 

We present a derivation of the orbital and spin sum rules for magnetic circular dichroic spec- 
tra measured by electron energy loss spectroscopy in a transmission electron microscope. These 
sum rules are obtained from the differential cross section calculated for symmetric positions in the 
diffraction pattern. Orbital and spin magnetic moments are expressed explicitly in terms of ex- 
perimental spectra and dynamical diffraction coefficients. We estimate the ratio of spin to orbital 
magnetic moments and discuss first experimental results for the Fe £2,3 edge. 
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Electron energy loss spectroscopy (EELS) in a transmission electron microscope (TEM) 
gives access at high energy losses to the density of unoccupied valence states with a sub- 
nanometer spatial resolution^. The possibility of using EELS to measure an energy loss 
magnetic chiral dichroism (EMCD) spectrum analogous to the X-ray magnetic circular di- 
choism (XMCD) signal obtained with synchrotron radiation^ has been suggested in 2003^ 
and demonstrated recently^. The principles of an EMCD experiment are the following: after 
suitable orientation of the sample, the incident electron beam is first elastically diffracted by 
the crystal. Each diffracted beam is then inelastically scattered. The total inelastic signal 
can be written as the sum of two kinds of contributions: the first one is due to each single 
diffracted beam and can be written in terms of the dynamic form factors (DFF) S (q, q, E), 
where q is the momentum transfer which depends on the diffracted beam and on the loca- 
tion of the spectrometer aperture in the diffraction pattern, and E is the energy loss. The 
second contribution involves all the possible pairs of diffracted beams and is described by 
the mixed dynamic form factors (MDFF) S (q, q', E)£. 

Inelastic scattering events are due to the coulomb interaction between the electrons of the 
probe and the sample electrons. In quantum electrodynamics, this interaction is described 
in terms of a virtual photon exchanged between the two electrons. The virtual photon 
associated to one of the DFF is linearly polarized in the direction of the transfer momentum 
vector (Lorentz Gauge). The polarization of the virtual photon associated to one of the 
MDFF is more complicated for any couple of transfer momentum vectors q and q'. It 
becomes right or left circularly polarized when the phase difference between the diffracted 
beams is ir/2, and when the momentum transfer vectors q and q' are orthogonal with 
identical modulus. This is the case for the vectors {q^q'i} and {q 2 ,q 2 } shown in Fig. QJi 
for which (q 1; q' x ) = n/2 and (q 2 ,q' 2 ) = — 7r/2. The EMCD signal is in this case obtained 
by subtracting the spectra measured at the two positions posl and pos2 shown in Fig. QJi. 
An accurate description of the EMCD spectra is not trivial, firstly because all the pairs 
of diffracted beams must be considered together, secondly because the propagation of the 
diffracted beams must be described within the fast electron dynamical diffraction theory, 
the incident and scattered electron beams behaving like Bloch waves inside the crystal. 

Magnetic circular dichroism has been measured in a TEM on the L 23 edges of 3d magnetic 
metals^. The most recent experimental papers describe the configurations which give the 
highest dichroic signal as well as dichroic/noise ratio. Several configurations have been 
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tested to reach this aim, like using convergent instead of parallel incident beam to increase 
the total current (LACDIF configuration^^) , or choosing the sample orientation and 
searching for the positions in the diffraction pattern which enhance the dichroic signal. 
This experimental investigation has been done by moving the diffraction pattern over the 
spectrometer aperture^, or with the energy spectrum imaging technique (ESI), which consists 
in recording the whole diffraction pattern for successive energy windows of typically 1 eV 
running over the L 2 3 edges^. The EMCD signal has also been calculated theoretically. 
These calculations, which are based on the first principles determination of the fast electron 
Bloch wave functions^ and transition matrix elements^ 1 ^ have been very helpful to find the 
experimental conditions and sample characteristics which give the highest dichroic signal. 

Up to now, EMCD experiments have only been analysed quantitatively in terms of dichro- 
ism. A quantitative interpretation of the spectra requires the determination of new sum 
rules which take into account the dynamical diffraction effects. In this paper, we present 
the analytical derivation of the orbital and spin sum rules for magnetic chiral dichroic spec- 
tra measured by EELS, and we discuss to which extent these sum rules can be applied to 
experimental results. 

The EELS spectra measured in a TEM can be described by the differential cross section 
d d ^Q n for scattering of a fast probe electron with the energy loss E and scattering angle f2. 
When the spectrometer aperture is located at a given position in the {x, y} plane of the 
diffraction pattern, the differential cross section for the core electron excitation edges can 
be written as^ 



o2 A det 

E% a S( q ,q,B) + E £ 2Re 



dEdn ^ <]* q qvq 



J^det 



q 2 q' 2 



(1) 



where the mixed dynamic form factors of one atom are given by 



S (q, q', = exp (-iq • r) |/) </| exp (iq' • r) \i) 5 (E - E f + E t ) . (2) 

\i) and |/) are the initial core states and the final unoccupied valence states with energies 
Ei and Ef.q = OS — g + q z e z and q' = OS — g' + q z e z are the momentum transfer vectors 
which depend on the vector OS connecting the transmitted beam and the spectrometer 
aperture in the diffraction pattern, on the reciprocal lattice vectors g and g', and on the 
momentum q z < which is transferred in the incident beam direction (Oz). The first and 
second terms in the right hand side of Eq. (CQ) describe respectively the contributions from 



the DFF and MDFF. The double sum over q and q' implies that the pairs of Bragg spots 
are not counted twice. The coefficients ^4q.q/ contain all the information on the Bloch wave 
eigenvectors and eigenvalues which is needed to describe the fast incident and scattered 
electron beams within the framework of the dynamical diffraction theory. They can be 
calculated as shown recently^. They depend on the momentum transfers q and q', on 
the location of the spectrometer aperture, on the atomic structure of the crystal, on the 
thickness and orientation of the sample, and on the location of the ionized atom inside the 
sample. For more complex unit cells, Eq. ([TJ) needs to be generalized by summing over the 
different atom species. The operators r + = x + iy, r_ = x — iy, and r = z can be used to 
express the mixed dynamic form factors within the electric dipole approximation as 

, QxQ'x + QyQy f s , 2 , - Qxq'y ~ q y q' x , , ( 

S{q,q,E) = u - [fi + + + q z fi + i y — - (3) 

where 

H + = J2\(i\r + \f)\ 2 8(E-E f + E i ), (4) 
hf 

^ = J2\(i\r-\f)\ 2 5(E-E f + E t ), (5) 

hf 

and 

Vo = Y,\( l \ r o\f)\ 2 HE-E f + E l ), (6) 

i,f 

z being the quantization axis. In the following, we have considered a four-fold diffraction 
pattern with distance g between Bragg spots and the two spectrometer aperture positions 
posl and pos2 which are indicated in Fig. CDb. Such a diffraction pattern can be observed 
with bcc Fe or fee Ni crystals oriented in the (100) zone axis. The reciprocal lattice vectors 
are given by g = nge x + mge y (n and m being the integers associated to each Bragg spot), 
and the momentum transfer vectors are written as q = g (5 — n) e x + g (e — m) e y + q z e z for 
position 1, and by q = g (5 — n) e x — g (e + m) e y + q z e z for position 2 (5 and e > being 
real). The difference and the sum between the EELS signals measured at the two symmetric 
positions of the spectrometer aperture described above are given by 
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\ dEm ) !' an< ^ ^he P a i rs °f Bragg spots (n, m) and (n', m') are 



not counted twice in the double sum. In the perfect zone axis configuration, this equation 
can be simplified using 

APOsl _ Apos2 /o\ 

n,m\n' ,m' n,— m;n',— m! ' \ I 



n,—m\n',—m' 



These equations remain valid in the systematic row configuration which is reached by tilting 
the sample around the (Ox) axis. This tilt modifies the value of all the coefficients A n>m . n i >m > . 
In particular, contribution from the Bragg spots which are not located on the diffraction row 
can be neglected and A nmn i tm i ~ if m ^ and/or m! ^ 0. The two beam case is obtained 
after a second tilt of the sample around the (Oy) axis. This tilt changes again the value of 
the coefficients A n m;n / m / which become small except if m = 0, m! — 0, and small n, n'. In 
practice, this facilitates the numerical calculation of the essential Bloch wave coefficients. 
Eqs. flZJ) and (jSJ) show that (cr 2 — <7i) is proportional to — To express (er 2 + 0"i) in a 
form which can further be used to derive the EMCD spin and orbital sum rules, we use the 



two additional approximations + ji- 



(0H 



jxJ) and no ps ~ + /x + ji- 



Thanks to these approximations, (a 2 + cri) becomes proportional to (fi + + /io + A*-)- The 
spin and orbital sum rules for an EMCD experiment can then be derived, using the sum 
rules which have been obtained by B. T. Thole et al. and P. Carra et al. to analyze XMCD 
spectra^^. The new EMCD sum rules can be written as 



and 
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J L3+L2 (a 2 + a 1 )dE '2 N h 
where (S z ) /Nh, (L z ) /Nh and (T z ) /Nh are respectively the ground state expectation values 
of spin momentum, orbital momentum, and magnetic dipole operators per hole in the d 



bands. The coefficient K contains all the information related to the dynamical effects. It 
can be expressed as 
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can be calculated for a very well defined geometry. It will depend on the excitation 
error of the incident beam, the specimen thickness, the detector position and aperture size. 
Moreover, in the experiment one never can achieve a perfectly parallel beam. Convergence 
and partial coherence of the electron source make the precise calculation of K untenable for 
the time being. Still, Eqs. j9]) and (TTOl) can be used to obtain 

f La (*2-vi)dE-2f L2 {a 2 -a l )dE = 4 (S z ) + 14 (T z ) 

Jl 3+ lA*2 ~ *i) dE 3(L.) 1 1 

free from any dynamical coefficient, sample orientation and thickness. Eqs. ((9j), ( fTOl) and ( TTTT ) 
apply to a single absorbing atom of the sample. The extension of the foregoing derivation 
from a four-fold symmetric diffraction pattern to the general case is straightforward. 

We now briefly describe the experimental applicability of the EMCD sum rules. Experi- 
ments were performed using the SACTEM Toulouse, a TECNAI F20 (FEI) equipped with 
a spherical aberration corrector (CEOS), an Imaging Filter (Gatan Tridiem) and a 2k*2k 
Camera (Gatan). An iron sample was used as a test sample. By combining the techniques of 
tripod polishing and ion milling, we prepared a large flat area which was electron transpar- 
ent. The magnetisation of the iron film is saturated in the (Oz) direction by the field of the 
objective lens pole piece. The sample was oriented in (110) two beam configuration and the 
electron diffraction pattern was recorded using the ESI technique performed with a 1 eV slit 
in an energy range of [645 eV, 745 eV] for a total of 30 min exposure time^. The diffraction 
pattern is taken using the LACDIF configuration^^ with a 7.8 mrad convergence angle 
which strongly increases the EMCD intensity and the signal/noise ratio. Post process correc- 
tions of isochromaticity and drift detected on the ESI data cube were applied using a home 
made software written in the scripting language of Digital Micrograph (Gatan). Finally, 
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EELS spectra are extracted using the ESI data cube, for the two positions OS = |e x ± ~e y 
located on the Thales circle which passes by the transmitted beam and the Bragg spot. Two 
circular apertures of semi-angle a = 4.2 mrad were used in the numerical integration, and 
the recorded spectra are shown in Fig. Eh- 

The difference between the two spectra gives the dichroic signal which is represented in 
Fig. 03. Our spectra have not been processed for removal of the background due to the 2p- 
state to continuum states transitions, because the aim of this paper is to demonstrate the 
feasability of the method. A quantitative analysis of spin/orbital moments would necessitate 
a better signal/noise ratio as well as more involved data treatment. We have applied Eq. (fT2i) 
to our experimental results, integrating the EM CD spectrum in the energy windows [705 eV, 
715 eV] for the L 3 edge and [719 eV, 729 eV] for the L 2 edge. Neglecting the contribution of 
the magnetic dipole operator, this measurement has given (L z ) / (S z ) = 0.18 ± 0.05. This 
result is higher but with the same order of magnitude than the values 0.124^, 0.088^, 
0.133^ and 0.086^ which have been obtained from neutron scattering data, gyromagnetic 
ratio or XMCD spectra. 

This comparison shows that EMCD is now on the way to giving quantitative magnetic 
information. Experiments do nevertheless deserve improvements, optimizing the angular 
and energy windows for integration in order to increase the still poor signal/noise ratio. 
Small background matching problems can also occur between the L 3 and L 2 edges. This can 
be seen near 715 eV where the dichroic signal does not perfectly vanish. These background 
problems are due to the fact that the non dichroic part of the signal is not perfectly the 
same at the two symmetric detector positions in the two beam case. In this case, Eq. ((H)) 
does not exactly describe the experimental configuration. This problem may be minimized 
by working at a higher voltage, in order to decrease the curvature of the Ewald sphere, or 
by looking for more symmetric experimental conditions for which Eq. JH]) holds perfectly. 

We have derived a set of sum rules for EMCD spectra which can be used to obtain orbital 
and spin moments of magnetic samples. Also dynamical diffraction effects of the electron 
beam in the specimen influence the dichroic spectra in a complicated way, the (L z ) / (S z ) 
ratio can be extracted straightforwardly when the scattering conditions are properly cho- 
sen. The main advantage of using EELS instead of X-ray absorption for this quantitative 
analysis comes from the subnanometer probe size which can be reached in a TEM. This 
opens exciting perspectives for the local magnetic analysis of nanomaterials and nanode- 
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vices like magnetic tunnel junctions for spintronics applications or magnetic nanoparticles 
with enhanced anisotropy and magnetisation. 
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FIG. 1: Diffraction pattern for an EMCD experiment, (a): experimental diffraction pattern for an 
iron sample oriented in (110) two beam configuration. The transfer momentum vectors for the two 
symmetrical positions posl and pos2 are represented by arrows, (b): Four fold diffraction pattern 
which has been used to express the differential cross section. The two different positions of the 
spectrometer aperture (OS = 5ge x ± eg^y, with 5 and e real numbers) which have been considered 
are indicated by open circles. The Bragg spots are represented by filled circles. A pair of integers 
(n,m) is associated to each Bragg spot, as shown for four of them. 
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FIG. 2: (a): EELS signal measured at the two symmetric positions posl and pos2 in the diffraction 
pattern of an iron sample oriented in the (110) two beam configuration; (b): corresponding dichroic 
signal 
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